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METHOD AND PROGRAMMABLE APPARATUS FOR QUANTUM 

COMPUTING 

5 

CROSS-REFERENCE TO RELATED APPLICATION 

This application claims the benefit of a provisional patent application no. 
60/447,566, filed on February 14, 2003. 

10 TECHNICAL FIELD 

The present invention is related to a method and programmable apparatus for 
quantum computing, and, in particular, a general method and programmable apparatus 
that uses quantum state propagation for quantum computing. 

1 5 BACKGROUND OF THE INVENTION 

Computer technology has evolved at an amazing rate in the past few decades. 
Further increases in the density of semiconductor integrated circuits may soon end 
due to the practical difRculties in patterning features with dimensions smaller than the 
wave length of ultraviolet light and soft x-rays. Ultimately, features cannot possibly 

20 be smaller than the length scale of single atoms and molecules. 

At the molecular level, quantum effects may provided huge advantages as well 
as problems not encountered at larger dimensions. Quantum computing has recently 
emerged as a possible approach to developing smaller computing devices that operate 
on the atomic level. In general, a traditional computer encodes information in a series 

25 of bits for computation that are manipulated via Boolean logic, which is based on 
binary mathematics. For example, a 32-bit digital computer generally manipulates 
one or several 32-bit operands during execution of a single instruction. A single 
instruction generally produces one or several 32-bit result values, each representing 
one of i?^ different possible result values. Thus, a single instruction maps one or 

30 several operands to one or several results, each within a range of 2^^ values. On the 
other hand, a typical quantum computer utilizes the +i4 and -14 nuclear spin states for 
isotopes such as ^^C and ^^F to represent the Boolean logic one and zero. In addition, 
the nuclear spin is also a quantum mechanical object that can exist in a superposition 
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of +14 and •V^ nuclear spin states. In theory, a typical 32-quantum-bit quantum 
computer can produce any combination of ip- values (i.e., 2^^! combintations) as a 
result of executing an instruction, providing for massive parallelism. However, in 
order to take advantage of the inherent Vi-spin state of isotaopes such as ^^C and *^F, 
5 the processor must be operated at temperatures near absolute zero. 

Computer manufacturers, engineers, and physicists have recognized the need 
for quantum computers that can handle the increasing demand for speed and 
computational volume and that are capable of operating at room temperature. 

10 SUMMARY OF THE INVENTION 

One embodiment of the present invention is directed towards a method and 
apparatus for quantum computing. In a described embodiment, data, a program, and 
uninitialized variables, or program output is provided. A compiler program converts 
the source code into an assembler code of computable functions. The computable 

15 functions are encoded by an encoding function. A continualization method is 
employed to determine a first-order, time-dependent, differential equation from the 
encoded functions. The variational principle is used to construct a Lagrangian whose 
minimum geodesic is the solution for the first-order, time-dependent, differential 
equation. The Lagrangian is converted into the quantum, canonical, Hamiltonian 

20 operator which is realized as an excitation field produced by an excitation generator. 
A control and scheduling system is employed to moderate the intervals between 
discrete applications of the excitation field produced by the excitation generator. The 
excitation field is repeatedly applied to a quantum processor consisting of a lattice of 
polymer nodes. The light radiation emitted by the quantum processor lattice of 

25 polymer nodes is converted by a transducer into an intensity-versus-vibrational- 
frequency spectrum. A running average of the intensity-versus-vibrational-frequency 
spectrum is maintained in the coherent memory. The average spectrum intensity 
values are used as coefficients in a polynomial approximating the encoding function. 
The output consists of the original data and program, and the substantiated variables, 

30 or program output. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 shows a control-flow diagram for the highest-level view of the 
method of one of many possible embodiments of the present invention. 

Figure 2 shows a control-flow diagram for a mid-level view of the method and 
5 apparatus that represents one of many possible embodiments of the present invention. 

Figure 3 shows a diagram of the apparatus comprising the "functional 
architecture" in Figure 2 that represents one of many possible embodiments of the 
present invention. 

Figure 4 is a control-flow diagram for the routine "functional architecture." 
10 Figures 5A-B illustrate a simple example input function. 

Figures 6A-B illustrate the application of the encoding in equations (3) - (6) to 
the simple example function. 

Figure 7A is an illustration of continualization of the encoding function F in 
Figures 3A-B via the step function in equation (8). 

15 Figure 7B is an illustration of a hypothetical interpolating function ^(y) . 

Figures 8A-F illustrate the continualization method for a hypothetical 
sequence in three-dimensions. 

Figure 9 is an illustration of a hypothetical geodesic solution x{f) of equation 
(23) on the three-dimensional carrier manifold. 
20 Figures lOA-C illustrate the Morse potential and spectrum for a hypothetical 

diatomic molecule. 

Figure 11 is an illustration of a quantum processor having a two-dimensional 
array of nodes representing one of many possible embodiments of the present 
invention. 

25 Figure 12 shows a simple cubic lattice of equally spaced nodes, reflective 

plates, and insulating boundaries representing one of many possible embodiments of 
the present invention. 

Figure 13 shows realization of the excitation field Hf generated by the 
excitation generator. 

30 Figure 14 is a diagram of the two-level automaton for a single node. 

Figure 15 illustrates the relationship between the chattering combination in 
equation (51) and the semi-open time intervals in equation (52). 
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Figure 16 illustrates chattering for a single hypothetical node using three pure 
states in the chattering combination (51). 

Figure 17 shows examples of lattice node interactions. 

Figure 18 illustrates the interactions characterized by the interaction 
5 Hamiltonian Hi. 

Figure 19 is an illustration of forward and backward propagating wave trains. 

Figures 20A-C illustrate a single hypothetical run of the excitation field 
Hamiltonian HfXxsQd to obtain a hypothetical spectrum from the quantum processor. 

Figure 21 shows example hypothetical spectra and the average spectrum after 
10 n hypothetical runs of flie excitation field Hf 

Figure 22 is a control-flow diagram depicting determination of the transducer 
read out period that represents one of many possible embodiments of the present 
invention. 

1 5 DETAILED DESCRIPTION OF THE INVENTION 

The present invention relates to a method and programmable apparatus termed 
flie Quantum Function Evaluator that uses quantum state propagation for 
computational processing. Figure 1 shows a control-flow diagram for the highest- 
level view of the method of one of many possible embodiments of the present 

20 invention. In step 102, the input consisting of data 104, a program 106, and 
uninitialized variables 108, or program output, is provided. The data 104 may be in 
the form of a separate file, input by an operator at the keyboard of computer, or 
actually incorporated into the text of the program 106. The program 106 may be a 
source code written in a high level programming language such as C, C-H-, 

25 FORTRAN, etc., or a lower level language such as assembler, or high level language 
such as scripting languages, shell scripts, or higher level mathematical programming 
languages. The uninitialized variables 108 represent program output, and may 
include data values transferred to any of many different output devices and types, 
including display terminals, electronic communications media, files or database 

30 management system data stored on mass storage devices, and many other types of 
output. Edge 1 10 represents executing one embodiment of the present invention that 
in a single iteration converts the input in step 102 into the output in step 112. In step 
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1 12, output includes the same data 104 and program 106. The uninitialized variables 
108, however, are now instantiated variables 114 determined by execution of the 
method and apparatus of the present invention. Note that if the data 104 in step 102 is 
changed, then a separate execution of steps 102 and 112 is needed to determine the 
5 instantiated variables 114. 

Figure 2 shows a control-flow diagram for a mid-level view of the method and 
apparatus that represents one of many possible embodiments of the present invention. 
Step 202 is the data 104, program 106, and uninitialized variables 108 described in 
step 102 in Figure 1. Steps 204 through step 216 provide a more detailed description 

10 of the method represented by edge 1 10 in Figure 1. In step 204, a compiler program 
converts the source code and data in step 202 into an assembly code representing a set 
of discrete valued computable functions. In step 206, the assembly code is converted 
into a discrete encoding of partial recursive functions. In step 208, a continualization 
method is employed to generate a continuous, first-order, time-dependent, differential 

15 equation from the encoded function determined in step 206. In step 210, the 
variational principle is employed to construct a Lagrangian whose geodesic is the 
solution for the first-order, time dependent, differential equation determined in step 
208. In step 212, the Lagrangian is converted into a quantum, canonical, Hamiltonian 
operator. In step 214, the quantum, canonical, Hamiltonian operator is physically 

20 instantiated by producing an excitation field described by the quantum, canonical, 
Hamiltonian operator. In step 216, the routine "functional architecture," representing 
one of many possible embodiments of the present invention, is employed to determine 
the coefficients of a polynomial approximation to the encoding function in step 206. 
After the polynomial approximation is determined in step 216, steps 206 through 202 

25 are performed in reverse order. In step 204, the polynomial approximation of the 
encoding function can be used to determine the computable functions in step 204, 
which leads to step 202, where the instantiated variables 114 in Figure 1 are given. 

Figure 3 shows a diagram of the apparatus comprising the "functional 
architecture" 216 in Figure 2 that represents one of many possible embodiments of the 

30 present invention. Located in the center of Figure 3 is a control and scheduling 
system 301 that controls and regulates the execution of a function input 302, an 
excitation generator 303, a quantum processor 304, a transducer 305, and a coherent 
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memory 306, as indicated by edges 307 - 31 1. Edges 312-317 indicate tfie direction 
in which information flows between the elements of the functional architecture. The 
function input 302 is employed to carry out steps 202 through 214 in Figure 2 to 
determine the excitation field and pass the Hamiltonian characterizing the excitation 
5 field to an excitation generator 303. The excitation generator 304 realizes the 
excitation field which is used to induce an emission of light radiation from the 
polymer molecule nodes of the quantum processor 304. The excitation field 
generated by the excitation generator 303 may be an electric field, light, radio waves, 
or any other means that can be used to generate vibrational excitation of polymer 

10 molecules. The transducer 305 converts the emitted light radiation from the quantum 
processor 304 into a coherent spectrum of intensities and corresponding frequencies 
that are passed to, and stored by, the coherent memory 306. The transducer 305 is also 
responsible for clearing the light radiation detector after each reading of the emitted 
light radiation from the quantum processor 304, as indicated by edge 318. The 

15 coherent memory 306 maintains a running average of the intensities and frequencies 
received by the transducer 305. The control and scheduling system 310 repeatedly 
executes the excitation generator 303, the quantum processor 304, and the transducer 
305 until the average spectrum intensities converge. Upon convergence of the 
spectrum, the average spectrum intensities are passed to read out 319 and 

20 subsequently used as the coefficients of the polynomial that approximates the 
encoding function in step 206 in Figure 2. 

Figure 4 is a control-flow diagram for the routine "functional architecture." In 
step 402, the function input determines the excitation field Hamiltonian according to 
steps 202 through 214, described above in relation to Figure 2. In step 404, a do-loop 

25 representing the control and scheduling system 310, described above in relation to 
Figure 3, executes steps 406, 408, 410, and 412. In step 406, the excitation generator 
realizes the excitation field Hamiltonian by bringing into physical existence an actual 
excitation field. In step 408, the quantum processor is impinged by the excitation 
field causing the polymer nodes to emit light radiation. In step 410, the emitted light 

30 radiation is converted into an intensity-versus-light-wave-frequency spectrum and a 
running average is maintained. In step 412, if the running average converges, then in 
step 414, the spectrum of intensities are used as coefficients of the approximating 
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polynomial to the encoding function in step 206 in Figure 2. In step 412, if the 
running average has not converged, then steps 406, 408, 410, and 412 are repeated. 

I. Quantum Computing 

A. Input Function and Encoding 

The computable functions generated in step 104 above are discrete valued 
functions referred to as "input functions." The input functions may be of the form: 



10 fiS^X'-xS^^D 

(1) 

where A^is the dimension of the domain 5, x--x5^, and the sets 5/ and D are finite 
subsets of the natural numbers given by: 

15 

5,={0,...,iV,.} 
D = {0,...,A^,} 

where Ni and Nd are the largest natural numbers in the sets 5, and D, respectively. The 
20 function /in equation (1) maps an A/-tuple denoted (xoy . . xn-i) in the domain 
S^X'"xS^, where Xi.j is an element in the set Si, to a single natural number in the 
range D. The AT-tuples (xq, . . x/^/.j) are referred to as "points" in the domain 
SiX'-xS^. 

Figures 5A-B illustrate a simple example input function having the form of 
25 equation (1) given by: 

f:S^xS2->D 



where the domain 5, is the set of all possible combination of points (xq, x/) such 
30 ihatxo andjc/ are respectively elements in the set: 
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5, ={0,1,2} 

and 

^2 ={0,1,2} 

The range of the example function /is given by the set: 

Z) = {0,1,2} 



Figure 5A shows a three column table of domain and range values of example 
function/ The ;co-column 501 and jcy-column 502 contain the domain elements of 
5, xSj, and the range values are located in /column 503. Figure 5B is a three- 
dimensional, Cartesian-coordinate depiction of the points (xq, xj,J) given in the table 

15 in Figure 5 A. The collection of points (xq, xjyJ)aTe plotted with respect to the jc^raxis 
505, the x/-axis 506, and the /axis 507, For example, the point 508 in Figure 5B is 
the graphical representation of the point in row 504 in Figure 5 A. 

It should be noted that, although two and three-dimensional illustrations, such 
as Figures 5A-B are provided in the following discussion, they are provided only for 

20 illustration purposes. In real-world applications, there may be many thousands, 
hundreds of thousands, or millions of variables, leading to extremely high 
dimensional problem domains. In general, these problem domains may be considered 
to be hyper-dimensional volumes, manifolds, or carrier manifolds. It is not possible 
to illustrate such problem domains, but techniques used to address hyper-dimensional 

25 domains may be analogized to three-dimensional illustrations. 

The input function may also be a solution of an iterative process of the form: 

yn^u=f^{ynA^'"^ynjt) (2) 



30 where i = 1 , . . A, and each 
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where 5= {0, 1, . . .,Ns}. 

The encoding of computable functions as described above in step 206 of 
5 Figure 2 can be accompHshed using the following function: 



10 



20 



F:[0,l,...,/7^-»]-^[0,l,...,p] (3) 



where 



X - 



= Z^. (4) 
and pis a natural number defined by: 
15 p = m^{{N,\i = h...,N},N,}+l (5) 



and 



^( .^UU,...^.-.) ifdefined 
^ ^ lO otherwise ^ ^ 



The encoding functions (3) - (6) map the domain values in iSj x- -XiS^y^ to a subset of 
points of the real numbers such that the values of F coincide with the values of/in the 
range of / Furthermore, the encoding functions (3) - (6) map x to 0 for any x that 
corresponds to an iV-tuple not in the domain of f. 
25 Figures 6A-B illustrate application of the encoding in equations (3) - (6) to the 

simple example function. The value for the natural number p given by equation (5) is 
given by: 
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/7 = inax{{2,2},2}+l = 3 

The encoding function in equation (4) maps the points {xo, xj) in the domain S^ x ^2 to 
a subset of the real numbers according to the expression: 

5 

1 

(7) 

In Figure 6A, the jc-column 601 are the jc-values calculated according to equation (7), 
10 and F-colximn 602 are the corresponding encoded function values F, In Figure 6B, 

the collection of points (x, F) displayed in Figure 6A are plotted with respect to the x- 

axis 603 and the F-axis 604. For example, the point (1, 2, 2) 508 in Figures 5B is 

mapped to the point (5, 2) 605 in Figure 6B. 

The example input functions given in equations (1) and (2) are functions of a 
15 larger class of discrete-valued functions generated by the compiler program in step 

204 in Figure 2. The class of functions generated in step 204 is referred to as the 

"Class r of computable functions and are given as follows: 

(i) Class / contains all the functions F encodable from a finite function/ as in 
20 equations (3) - (6). 

(ii) If 

{/i - » A : ^. X • • ' X 5^.1 , ^ finite, 5, = {l, A^, | A^, finite}} 

25 

is a subset of Class /, so is the direct sum: 

y;©..-©/^:5,x...x5^->5;;^^, 



30 



(iii) Class / contains the projection functions: 
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where p, (n,,...,w,,...,/i^) = ai, 

5 (iv) If 

{f,,^.-.f^\frS,X'-xS^^S^,,,Kfmite,S, = {^^^^^^ finite}} 
is a subset of Class / and 

10 



is in Class /, so is the composition: 
15 g(^...,/^):5,x...xS^->S^,, 
(v) If for each neN , 

g:NxS^x- -xSg -^N 

20 

is in Class / so is the following: 

min {«, g{n,n^, .... n^) = 0] 

25 (vi) If 

g-.NxS,, -^S^ 



is in Class /, then the family of functions given by: 

30 
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is in Class /. 

(vii) Any fiinction constructed by the finite application of (i)-(vi) is in /. 
5 The elements of Class / are either an encoding or can be transformed into an 

encoding function via a finite number of steps (ii) - (vii) above. Clauses (ii) - (iv) are 
included because these clauses are constructed by the compiler program found in 
many high-level programming languages. Although the Class / of computable 
functions includes only completely specified functions, in alternate embodiments the 
10 method of the present invention can be modified to handle the larger class of discrete 
partially specified functions and indicator functions. 

The central objective behind the method and apparatus of the present 
invention is to find effective and fast means of computing the discrete Class / 
functions via quantum approximations to the encoding functions. 

15 

B. Continualization 

The method of continualization is described in detail in patent application no. 
10/693,729 and is incorporated by reference. The example function provided in 
Figures S and 6 can be continualized using a step function O defined by: 

20 

O:[0,...,/7^-^)^[0,...,p^-^) (8) 
0{y) = F{x) for each:,;, e [;c,x + l),:,; g {0,...,p^-*} 

Figure 7 A illustrates continualizing the encoding function F shown in Figures 6A-B 
25 via the step function in equation (8). Continualization increases the size of the 
domain by including the real numbers between the discrete domain values x and 
assigns to each y a function value 0{y) . Figure 7A shows the graph of step-function 
<!>(;;) coincide respect to the x-axis 701 and the <I>-axis 702. The function values 
F(x) equal to ^{y) where the discrete x values equal the discrete domain y values. 
30 For example, the point (5, 2) 703 in Figure 7A is the point (5, 2) 605 in Figure 6B. 
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The continualization process is not limited to the step function given by 
equation (4). In alternate embodiments, other interpolation schemes can be used to 
construct a function 9{y) whose values coincide with F(x) such as linear 
interpolations, Adams-Boshworth interpolation, splines, Fade interpolation, etc. 
Figure 4B is an illustration of a hypothetical interpolating function ^{y)y where 
0(y) equals the encoding function F{x) for values of equal to the discrete 
domain values x. 

The interpolating function 0(y)m3y be expressed for hyper- 
dimension manifolds as a sequence of real vectors {jCq, jCi,...} defined by the 
following recursion relation given by: 

(9) 

where each vector is a A:-tuple in ^-dimensional set of real numbers given by: 

and sequence {a„] is a sequence of real numbers having the properties: 

lima =0 

n 

and h is function mapping from J?* to i?*. In other words, a new state-vector value 
x„+j is computed, in each iteration of a computation, from a current state-vector jc„ and 
a discrete, vector-valued function h that depends on x„. Figure 8A-F illustrate the 
continualization method for a hypothetical sequence in three-dimensions. Figure 8A 
shows the vectors jfn-2, ^^^n-i, Xn, and Xn+i plotted in with respect to the xy-axis 802, 
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the X2-axis 804, and the -rj-axis 806. Next, the sequence {x„ )^^^ can be approximated 
by the sequence {X„ having terms satisfying the recursion relation given by: 

5 (10) 

where b„ is a vector parameter. Figure 8B shows the sequence given by the 
approximating recursion relation in equation (10). The sequence of vectors {x„)^^^ 

shown in Figure 8A have been replaced by the sequence of vectors {X„)„^q in Figure 
10 8B. 

The time sequence can be defined as follows: 



(11) 

15 where =t.^^ 



Figure 8C is an illustration of the relationship between the time sequence [t^] and the 
real number sequence {a^} given in equations (9) and (10). Note that because each 
element a„ approaches zero as n approaches infinity, the real number sequence {a„} is 
20 the length of the decreasing time intervals -t^) on the time axis 808. For 
example, real number ay 810 is a longer interval than a„Sl2, 

Next, two time dependent vector functions X^{t) and g^{t,a,b) are defined 
as follows: 



25 



X''{t„) = X„ for all n^O 
(12) 
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(13) 

and 

{t„ ,a,b) = J^ a,b, for all « > 0 

/=0 

(14) 



(15) 



Equations (13) and (15) are linear functions with respect to time /. Figure 8D is an 
illustration of the functions defined by equation (13) and shows that segment {t) 
814 is a vector between the vectors Xn and Xn+i in R^, 
A function X^ can be defined by the following: 

(16) 

Figure 8E is an illustration of the equation (16) plotted as a function of time t. The 
function values are plotted as a step function to show that X^ (/) 816 is constant on 
the time interval (t^ , t^^^ ) . 

Summing the terms in equation (10) fi-om zero to /i - 1 gives the following 
expression: 

/=0 1=0 1=0 1=0 

(17) 
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After some algebra and applying the definition of the integral from integral calculus, 
equation 17 can be given as continuous, vector-valued function of time t as follows: 



X'^ (0 = X« (0) + \h(x' {s))ds + g'' {t,a,b) 



(18) 



The functions defined as: 



X"{t) = 



_{x\t+t„) iit>-t„ 



if t<-t. 



10 



can be used to rewrite equation (18) in the following form: 



X"{t) = 



X" (0)+ \h(x" {t„ +s))ds + g° {t,a,b) iff > -t„ 



ift <. -t. 



15 



(19) 



Equation (19) can be expressed equivalently as follows: 



X"{t) = 



X" (0)+ \h(X" {s))ds + g' it,a,b) + e" (t) ift > -t„ 



X„ 



if t^-i 



20 



where the sequence of functions e" (t) and g" {t,a,b) approach zero as n approaches 
infinity uniformly on finite time intervals -t„). The time dependent sequence 
{X" (f)} is bounded and equicontinuous on interval (-00, 00). Thus by the Arezla- 
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Ascoli Lemma, there is a convergent subsequence of the sequence of functions 



The sequence of functions {g"(t^a,b)} are also bounded and equicontinuous 

on interval (-00, 00). Again, applying by the Arezla- Ascoli Lemma, there exists a 
convergent subsequence of functions |^"(/,a,6)| whose limit is the vector function 

Thus the original discrete sequence given by equation (9) can be modeled 
by the first-order, time-dependent, differential equation given by: 



Figure 8F is an illustration of the differential equation (20) as a tangent vector 818 to 
the vector solution functions x(t) 820. 

C. Variational Model 

The variational model is used to formulate a Lagrangian for the first-order, 
time-dependent, differential equation (20) in step 210 in Figure 2. The procedure 
begins with the variational problem formulated in terms of the Lagrangian given by: 




X{t) = h{X{t)) 



where AT is a shorthand notation for the time derivative of X. 



x{t) = h{x{t))-g{t,a,b) 



(20) 



First, the Lagrangian is assumed to have to the following properties: 
(1) The Lagrangian is homogeneous of degree 1 in i: . In other words. 



L (x, yljc, t) = AL{x,x,t) for a real valued A > 0 . 
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25 



(2) The Lagrangian is positive definite in i . In other words, 



5 



dXjdx. 



Properties (1) and (2) allow for the definition of a metric ground form ds on 
the carrier manifold such that 

^^ = T^gij{x>xj)cix,dxj 

ij 

10 (21) 

A geodesic is the path that represents the shortest path between any two points when 
the path is restricted to a particular surface. The geodesies subject to the metric 
ground in equation (21) are the extremum of the corresponding parametric calculus of 
variations problem which is to find a curve on the manifold which minimizes the 

1 5 following integral : 

)L{x{t),x{t),t)dt 

to 

(22) 

20 satisfying the endpoint conditions. The geodesies are solutions to the Euler-Lagrange 
equation given by: 



d dL{x,xj) dL{x,xj) ^ 
ds dx. 5bc,. 



(23) 



Figure 9 is an illustration of a hypothietical geodesic solution x{t) 904 of equation 
(23) on the three-dimensional surface carrier manifold 902. The path followed by the 
geodesic solution j^t) 904 is the minimum path between the starting point x(td) 906, 
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which represents the data 104, program 106, and uninitialized variables 108 in step 
102 in Figure 1, and the ending point x(tj) 908 represents the data 104, program 106, 
and instantiated variables 1 14 in step 1 12 in Figure 1 . 

Second, the vector function jc(t) defines a trajectory that has an integral within 
5 an f of its minimum and is referred to as the ^-optimal control. 

Third, the Lagrangian L(x, u, t) is replaced by L*(x, u, t) which is the 
convexification of the Lagrangian L with respect to u. Then an existence theorem for 
the relaxed controls can be applied to L*(x, u, t) to ensure that the optimal control 
problem has a relaxed solution. The relaxed solution can be approximated to obtain 
10 an f-optimal control for the original problem. 

Fourth, for any pre-specified £ , the ^^ptimal control can be computed for the 
original problem and implement ^--optimal control as a finite-state, physically- 
realizable, control automaton. The actual control law issued by the finite-control 
automaton at state x \s 3. chattering control. The chattering is among approximations 
15 to some local extrema of L(x, u). 

The desired evolution of solution x[t) to the first-order, time-dependent, 
differential equation (20) is assumed to by define a path along a geodesic of the 
carrier manifold. In other words, there is a desired Lagrangian L{xyX,t) which 
induces a metric on the carrier manifold as described in equation (21) by the evolution 
20 of a function x{t) satisfying equation (20). The covariant partial derivatives on the 
carrier manifold are defined as follows: 



25 There exists a covariant partial derivative with respect to time t, however, time is an 
invariant, therefore the covariant partial derivative with respect to time t is just: 




and 




dj dt 
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Therefore, the transformation from local coordinates to coordinates on the carrier 
manifold is given by: 

where T/. are the Christoffel coefficient arising from the Levi-Civita connection 

associated with the metric in equation (21). The Christoffel coefficients may be 
computed explicitly via the formula: 

where (g'*^ ) is flie inverse of the metric ground form matrix given by: 



V ' ' ) 

15 

For convenience, the covariant partial derivatives The derivation begins by taking the 
derivative with respect to time t of both sides of equation (20) to give the following 
expression: 



20 



x = -r-x-g{t,a,b) 
a.x 



(24) 



where x = [xg] and x = [xg] are column vectors and 



25 



d.x 
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is a k-ak matrix. The Lagrangian Z, is a function mapping defined as follows: 

where TM is the tangent bundle and the Lagrangian L{^x,x,t) satisfies the Euler- 
Lagrange equation given by: 



d_ 
dt 



10 



(25) 



where — 
dt 



and — ^ are colimin vectors 



Expanding the time derivative in equation (2S) gives tiie following: 



15 



dL 
d,xd,x 



■x + 



dlL .. d 

£ ;c+— 

5,ia,i dt 



d^x 



= 0 



(26) 



The Lagrangian L{x,x,t) is assumed to be three times differentiable, and the mixed 
partial derivatives of L[x,x,t) are assumed to be continuous so that order of 
20 differentiation does not make any difference for mixed partial derivatives of 
L{x,x,t) of degree 2 or 3. Thus taking the partial derivative of equation (26) with 
respect to of the ith component gives the following: 

^ dlL dlL ^ dlL f^d^hj . .. , ^."1 

jd,x,d^x,d^Xj d^x,d^x, jd^x,d^x,d^Xj\y dx ') 

25 (27) 
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Equation (27) can be rewritten as: 



^ dlL ^ 



dt 



(28) 



dlL dh dlL dlL 



\^d^x,dxj yd^Xjd^x, d^x, d^x.d^x, d^x.d^x, 



^ = 0 



Thus, equation (28) can be rewritten in matrix form as follows: 



d_ 
dt 



f a2 



dlL 1^ dlL dlL 



\^d^xd^x) d^xd^x d^xd^x 



^d^xd^x^ 



dh 



d,x 



(29) 



Taking the transpose of both sides of equation (29) and using the partial differential 
property given by: 



SitL 



dlL 



d^xd^x) yd^xd^x^ 



gives the following: 



dQ 
dt 



where 



25 



Q = 



d^xd^x 



(30) 
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Equation (30) can be explicitly solved for by integrating along the path of the 
characteristic, second-order, partial-differential equation given by: 

5 

D. Quantization 

In this section the development turns to an alternative statement of the 
structure and formulation known as the "quantum, canonical, Hamiltonian operator" 
as first introduced in step 212 in Figure 2. The usefulness of the Hamiltonian 

10 viewpoint lies in providing a framework for theoretical extensions in many areas of 
physics and provides much of the language with which present day quantum 
mechanics is constructed. 

In tiie Lagrangian formulation given above, the system is characterized by the 
time-dependent, vector-valued function x{t) of k independent degrees of freedom 

1 5 and is a problem with k independent variables jcjt, and k time derivatives Xj^ . The 
Hamiltonian is fundamentally different from the Lagrangian, because the Hamiltonian 
is formulated in terms of the position coordinates and the associated momenta pk 
(Classical Mechanics Edition, Herbert Goldstein, Addison- Wesley Publishing Co., 
New York, 1980). The Lagrangian formulation is transformed into a Hamiltonian 

20 formulation by first taking the total time derivative of the Lagrangian L[x,xj) to 
get: 

dL ^ dL . ^ dL .. dL 
dt Ydx, * ^dx, * dt 

(31) 

25 

Next, the Euler-Lagrangian equation (25) is substituted into equation (31) to give the 
following equation: 
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dt 



(32) 

The quantity in the brackets is referred to as the "energy function" and is denoted by: 



(33) 



The energy function (33) can be reformulated in terms of the linear conjugate 
10 momenta pk which are associated with the coordinates Xk and are determined by the 
equation: 



dL 

Pk ^-^T 



(34) 



1 5 The conversion from the Lagrangian formulation to the Hamiltonian formulation is 
treated strictly as a mathematical problem by substituting equation (34) into equation 
(33). Thus, the energy function given by equation (33) can now be represented as a 
Hamiltonian function of the coordinates Xk and conjugate momenta pk given by: 



20 H{x,pj) = '£p,^x,-L{x,x,t) 

k 

(35) 

Equation (35) is written in this manner to stress the fact that the Hamiltonian is 
always considered as a function of the coordinates {x^pj) , whereas the Lagrangian 
is a function of the coordinates (jc,jc,/) . 
25 Using the classical mechanical expression for the Lagrangian given by: 
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L{x,x,t) = Y,^m,xl-V 

k 

where F is a potential energy function of x or both x and time t, equation (35) can be 
5 rewritten as a function of the canonical variables (jc, p) to give the classical 
mechanical Hamiltonian: 

k ^^k 

(36) 

10 

If the potential energy is firee of the time t, the system described by the Hamiltonian in 
equation (36) is said to be conservative, which means that no energy is being added to 
or removed from the system. 

The Hamiltonian expressed in equation (36) best characterizes dynamical 

15 systems that behave in a continuous manner. However, the Hamiltonian as 
formulated in equation (36) is not adequate for characterizing the quantum 
mechanical method and apparatus of the present invention, because the present 
invention reUes on the behavior of atoms and molecules, which exists in a certain 
number of stationary or quantum states. 

20 The Hamiltonian in equation (36) can be modified to characterize the quantum 

behavior exhibited by atoms and molecules by substituting the differential operators: 

dt 

25 and 

h d 
Pk^-^ 
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into equation (36) {Quantum Mechanics, Albert Messiah, Elsevier Science Publishers, 
The Netherlands, 1961). The resulting equation is the quantum, canonical, 
Hamiltonian operator given by: 



\ i dx dt^ 



= > = //i — 

) V dxl dt 



(37) 



that characterizes the behavior of a quantum mechanical system such as an atom or 
molecule. Equation (37) is called an "operator** because the kinetic energy and total 
10 energy have respectively been replaced by the differential operators given by: 



3' ^3 
— and — . 

dxl 



The Hamiltonian in equation (37) is a second-order, differential, wave 
15 equation. In general, the solutions to the Hamiltonian given in equation (37) are 
complex-valued, wave functions that depend on spatial coordinates, time, and integer 
parameters. The wave fiinction solutions of the Hamiltonian are given by: 

20 (38) 

where the integer parameter n is referred to as the "quantum number." In general, 
there are an infinite number of solutions referred to as "eigenfunctions," 

"eigenstates," or "states," that contain all the information that can be determined 
25 about the state of the system such as position, momentum, and energy. The set of 

solutions II given by equation (38) form a basis for a Hilbert space, and 

therefore, given any two eigenfunctions | ^„ (t)) in the Hilbert space : 
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(39) 



where * is the complex conjugate of ^, . 

5 The Hamiltonian given by equation (38) can also be rewritten as an eignevalue 

problem in the form: 



10 



15 



(40) 



where A„ is the energy of the system in the state |*P„(^)) and is referred to as an 

"eigenvalue." Prior to a measurement of the quantum mechanical system 
characterized by the Hamiltonian, the system is assimied to be represented by a linear 
combination of state functions |T„ (/)^ as follows: 



(41) 



where is the state function for the entire system and d are complex- valued 

20 coefficients. When a measurement of the system is made, the system is observed in 
only one of the eigenstates in equation (41), where the probability of 

observing the system in the state \ is the square of the coefficient, \c.f . 

The potential energy function V in equation (40) characterizes the particular 
physics of the system under observation such as the rotational and/or vibrational 
25 motion of a molecule, or the electronic states of an atom or molecule. For example, 
substituting a potential V characterizing the vibrational behavior of a diatomic 
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molecule, referred to as Vvit, into the Hamiltonian, gives the vibrational eigenstates 
(^)) accompanying vibrational energy eigenvalues ^. 

Figures lOA-C illustrate the Morse potential and spectrum for a hypothetical 
diatomic molecule. The vibrational potential of a diatomic molecule may be 
5 approximated by the Morse potential given by: 

f'w.=A(l-exp(-vr))' 

where De is the energy required to completely disassociate the diatomic molecule, k is 
10 the curvature at the potential minimum, and r is bond length 1002 between atoms 

1004 and 1006 shown in Figure lOA. The Morse potential characterizes the 

interaction between the atoms 1004 and 1006. 

Figure lOB is an illustration of the Morse potential and the quantized 

vibrational energy levels. The Hamiltonian representing the vibrational behavior of a 
1 5 hypothetical diatomic molecule is given by: 



where is the set of vibrational states of the diatomic molecule, and {X^ is the 

20 set of corresponding eigenvalues. The potential energy Kv,^ is plotted as a function of 
the inter-atomic bond distance r 1002 of atoms 1004 and 1006 in Figure lOA, where 
the energy is represented by the vertical axis 1008, and the inter-atomic bond distance 
r is represented by the r-axis 1009 in Figure lOB. The lowest vibrational state is 
with the accompanying energy eigenvalue ^ 1010 is referred to as the "ground 

25 state." 

In Figure lOB, the light sources 1012 and 1014 emit photons with energy 
given by: 



Docket No. 35021.005 



29 

where v. ^ is the light frequency, and 

The hypothetical diatomic molecule absorbs a photon from source 1012 having 
5 frequency vojy and the diatomic molecule experiences an increase in vibrational 
energy by an amount \^ y as indicated by edge 1016. The transition from the ground 
state to a higher energy state is referred to as "excitation " and the molecule is said to 
exist in an "excited state." The hypothetical diatomic molecule now has an energy 
eigenvalue Aj and the state of the molecule has transitioned from the groimd state 

10 \zo) to the excited state \zx) • Edge 1018 shows a second transition from the state 
l^'i) to the excited state \z2) resulting from absorption of a photon with frequency 
vj,2 emitted from light source 1014. In this hypothetical example, the excited states 
are generated using light of a particular frequency. It should be noted that atoms and 
molecules may also experience excitation to higher quantum states by application of 

IS an electrical field. 

Atoms and molecules also emit photons to the surroundings as they transition 
from an excited state to a lower energy state. The process of transitioning to a lower 
energy state is referred to as "relaxation" or "thermal relaxation." Edge 1020 shows 
the relaxation of the hypothetical diatomic molecule from the excited state Izz) with 

20 energy eigenvalue A2 to the state | Zi) with energy eigenvalue A/. The transition from 
the state Izj) to the state \zi) results in the emission of a photon with the frequency 
vij. Edge 1022 shows a transition from the state \zi) to the ground state \zo) and 
the emission of a photon with the frequency vo.i- 

The emitted photon can be observed by a detector to generate a spectrum of 

25 the different observed frequencies of light for an ensemble of diatomic molecules. 
Figure IOC is an illustration of the ensemble average of the spectmm of gas phase 
hypothetical diatomic molecules as described in relation to Figures lOA-B. The 
horizontal axis 1024 is the observed quantized vibrational frequencies for the 
ensemble and the vertical axis 1026 is the intensity associated with each vibrational 
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frequency. The spectrum of the hypothetical molecule reveals a tapering off in the 
intensity with higher vibrational frequencies. The intensity is a measure of the 
fraction of diatomic molecules in a particular vibrational state. Emission line 1028 
has the largest intensity So and can be used to approximate the fraction of hypothetical 
diatomic molecules in the state \x^' Therefore, according to equation (41) the state 
fiinction describing the ensemble of diatomic molecules can be given by: 



E. Realization 

Realization is the process of transforming the canonical, quantum, 
Hamiltonian operator in equation (40) into a physical operation of the hardware to 
solve for the ftmction ^{x) . The fimction ^{x) may be expanded in terms of an 
orthogonal basis set of /7-valued Chrestenton frinctions spanning the Hilbert space Yp 
given by: 



where |c,.|^ ^5. 



where the general form of a Chrestenton basis frmction is given by: 




with 
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and 

J-Z^Js P 

5=0 

5 

Therefore, the interpolating functions ^{x) , written in terms of the orthogonal basis 
set Bp is given by: 

*s[0,/.''-'] 

10 (42) 
where 

15 (43) 

The set of coefficients is referred to as the spectrum of the function in 
equation (42), where the number of terms in the spectrum { } ranges from about 10^ 
to 10^^. If the spectrum {>8^} is known, then the function ^{x) can be used to 
20 determine the encoding fiinction F{x) as described above in relation to steps 216 and 
206 in Figure 2. The encoding function F{x) can then be used to determine the 
original function X^) as described above in relation to steps 204 and 202 in Figure 2. 
The objective of the present invention is to determine a spectrum that approximates 
the spectrum {A} . 

25 The hardware of the system is characterized by the hardware Hamiltonian 

operator referred to as "//i?." The realization step consists of constructing the 
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excitation-field Hamiltonian Hf, referred to as so that the hardware, when 

interacting with the excitation field Hf, satisfies the following condition: 

\H^^H^-H\<e 

5 (44) 

where | ■ | denotes an operator norm, and £* is an a parameter chosen to satisfy the 
precision requirements for the computation. 

In general, the eigenvalues of the composite Hamiltonian Hq + i/^ approximate 
10 the spectrum {p^^ of the function given by equation (42). Let \k = 1,...,//} be the 
subset of the eigenvalues of Hq + Hf with the corresponding eigenstates 
{|T^(^))|A: = l,...,iv} so that if 

15 (45) 
then 

(46) 

20 and 

^\4>.{x)-\^,{x))\lx<e 

(47) 



25 where is a Chrestenton basis function. Therefore, the interpolating function <I>(;c) 
given in equation (42) can be approximated as follows: 
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k 

(48) 

Thus, a discrete function can be approximated by exciting the hardware 
5 system Hq^ according to the excitation-field Hamihonian Hf to obtain the set of 
coefficients for equation (48), The encoding function can then be constructed fi-om 
equation (48) which in turn can be used to determine the Class / functions, as 
described above in relation to reversing steps 206 - 202 in Figure 2. 

10 II. The Functional Architecture for Quantum Computing 

A. Function Input 

The function input 302 in Figure 3 symbolically determines the excitation- 
field Hamiltonian /(^via the steps 202 - 214 described above in relation to Figure 2. 
The process of symbolically determining the excitation-field Hamiltonian Hf is 
1 5 referred to as the "compilation process." 

B. Excitation Generator 

The excitation generator 303 in Figure 3 receives the symbolic excitation-field 
Hamiltonian Hf bom the function input 302 and converts the symbolic description of 
20 Hf into a physical implementation to be carried out by the quantum processor. The 
excitation generator addresses each node of the quantum processor according to the 
desired state behavior given by: 

25 (49) 

where | (/)^ is the state function of the yth node in the quantum processor. 

C. Quantum Processor 
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Figure 11 is an illustration of a quantum processor 304 in Figure 3 having a 
two-dimensional array of nodes representing one of many possible embodiments of 
the present invention. The nodes 1102 are chemical polymers linked to adjacent 
nodes in the lattice via quantum lateral bonds 904 and quantum forward bonds 1 106. 
5 The boundaries 1106 and 1108 provide an insulating barrier to prevent the electric 
field fi-om leaking out to the surrounding components of the system. The quantum 
processor also contains a excitation filed 1112 that implements the excitation-field Hf. 
Located at opposite ends of the quantum processor are reflective plates 1114 and 1116 
that reflect the excitation field back onto the lattice of nodes 1 102. 

10 The lattice of nodes is not limited to flie rectilinear nodal arrangement shown 

in Figure 11. In alternate embodiments, the nodes may be arranged to provide a 
higher density of nodes, as, for example, by offsetting the nodes in adjacent rows to 
produce a more closely packed arrangement of nodes. In alternate embodiments, the 
nodes may be arranged in any nimiber of different three-dimensional lattice structures, 

15 as, for example, the simple cubic, body-centered cubic, face centered cubic, primitive 
face centered cubic, simple hexagonal, hexagonal closest packed, rhombohedral etc. 
Figure 12 shows a simple cubic lattice of equally spaced nodes 1202, reflective plates 
1204 and 1206, and insulating boundaries 1208 and 1210 that represents one of many 
possible embodiments of the present invention. 

20 The Hamiltonian realizing the quantum processor prior to receiving the 

excitation-field Hf from the excitation generator is given by: 

H = H^ 

25 where Ho is the hardware Hamiltonian. Prior to application of the excitation field, the 
lattice nodes are in the ground vibrational state. After the quantum processor receives 
the excitation-field Hf fi-om the excitation generator, the quantum processor is 
characterized by the Hamiltonian given by: 



30 



H^H^+H^ 
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Figure 13 shows a hypotheitical excitation-field H/ applied to the hypothetical 
quantum processor shown in Figure 11. The excitation-field H/ originates at the 
excitation field 1302 and sweeps across the quantum processor in the direction 1304. 
The excitation field wave firont 1306 passes over nodes in their ground state exciting 
5 each node to a higher vibrational and/or electronic state. The nodes in the ground 
state are identified by the open circles 1308, and tiie nodes in an excited state are 
identified by tfie solid circles 1310. The rotational quantum states of the nodes are not 
accessed because each node is bound in the lattice and cannot fireely rotate. 

Figure 14 is a diagram of the two-level automaton for a single node 1402. The 

10 lower automaton 1404 provides a node's view of the quantum processor consisting of 
the transition selector 1406, the state transition 1408, and current state 1410. The 
upper automaton 1412 represents the quantum computing paradigm as described 
above by the fimction input 302 in Figure 3 and consists of the probability state 1414 
and the probability transition 1416. The input selector fimction 1418 is an interface 

15 fimction, carried out by the excitation generator 303 in Figure 3, that models the 
interaction of the node 1402 with the excitation field 1418, the neighboring states 
1420 of the other nodes in the lattice, probability states 1414, and the current state 
1410. 

20 1. The Lower Level Automaton 

In the lower automaton 1404, the state transition 1408 characterizes the 
programmable discrete spectrum of the node and is controlled by the mixing state 
probability density transition computed in the upper automaton 1412. The transition 
selector 1406 in the lower automaton 1404 executes the commands issued by the 

25 upper automaton 1412. Each node 1402 has a set of pure states given by: 

(50) 

30 where |^^^ (/)^ are the pure states of the node, J is the index of the yth node in the 
quantum processor, k is the quantum number, and Np is number of available quantum 
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states. Each node in the lattice begins in the ground state | {t)j and is excited to a 
higher states |^^^ (/)^ by the excitation-field H/, as described above in relation to 

Figure 1 1 . After the node has been excited to a higher state by the excitation-field H/, 
the node spontaneously undergoes relaxation to a lower energy state in S. After a 
5 period of time determined by the control and scheduling system 301 in Figure 3, the 
lattice of nodes are subsequently re-excited by another excitation-field H/ generated 
by excitation generator 303. The time duration between the excitation fields is 
referred to as the "update time" and is denoted by A. The update time A of the node is 
determined to be larger than 10 times the maximum relaxation time of any of the state 
10 transitions in the lower automaton. 

Each node in the lattice may transition between one or more states in S during 
the update time A and referred to as "chattering." The chattering effect is described 
as a combination of the pure states in S and is given as follows: 



l^.(-))= 



(51) 



where 



20 



(52) 



are the semi-open intervals of the update time interval [/, r + A) and give the duration 
that nodey spends in the pure state |^^^ (r)) . The time interval A^^ (t) is a fiinction 
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of the characteristic excitation or relaxation times associated with the state. See 
Fields and Nodes: Field Theory and Dispersion Relations^ by K. Nishijima, 
Benjamm, Inc., NY, 1969. 

Figure 15 illustrates the relationship between the chattering combination in 
5 equation (51) and the semi-open time intervals in equation (52). The pure states 
*® chattering combination (51) are arranged vertically in order of 



combination (51) and defined in equation (52) are displayed graphically below the 



1502 for the length of time A^.^ 1504 in the interval Ij^ (/) 1506. The sum of the time 
intervals A^. can be expressed as follows: 



increasing energy. The semi-open intervals Ij^{t) appearing in the chattering 



10 



pure states y¥ {tyj . Each semi-open interval Ij^ {t) is of length A^^ and represents 
the amount of time the yth node remains in the pure state j^^.^ {t^ during the update 
time interval t +A). For example, the 7th node occupies the pure state |^^. {t^ 



15 



k 



(53) 



The chattering coefficients can then be defined as: 



20 




A 



(54) 



which allows equation (53) to be written as: 



25 



k 



(55) 
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A chattering node may spontaneously transition to a lower state in S giving off 
radiation to the surroundings, or the node may absorb radiation emitted by other nodes 
in the system causing the node to transition to higher energy states in S. Figure 16 
5 illustrates chattering for a single hypothetical node using three pure states in the 
chattering combination (51) arranged in order of increasing along the energy-axis 
1602. First, assume the node begins in state [^^^^ (r)^ 1604 with energy Aj^ 1606. 

The node in the state j^^^^ (r)^ 1604 spontaneously transitions to the lower state 

I^A-i (^)) ^^^^ energy Aj^ ^ 1610 emitting energy hv.^ ^ j^ to the surroundings, 

10 as indicated by edge 1612. Second, the node absorbs energy hVj^^j^ from the 
surroundings, as indicated by edge 1614, the node transitions back to the state 
I^A (^)) with energy 164. Third, the node absorbs energy with frequency 
hVj^j^^^ from the lattice surroundings, and transitions to ttie state |^^.^ ^ (r)^ 1616 with 
energy Xj^^^ 1618, as indicated by edge 1620. 

15 For each update time interval [/,/ + A), the transition selector 1406 receives a 

command from the upper automaton 1412 in Figure 14, which is comprised of an 
ordered tuple of chattering coefficients jor^^ (0|A ^^p} • The Transition Selector 

1406 then computes the chattering combination of mixed states {t^ according to 
equation (51). 

20 

2. The Upper Level Automaton 

In general, there is not enough information to determine a specific state 
function that characterizes the lattice or any of the lattice nodes. Therefore, the best 
description of the computation is a probabilistic description. In the quantum 
25 formalism, the description is referred to as the "probability density description." See 
Probabilistic Model of Quantum Relation, V. F. Weisskopf and E. Wigner, Physik 63, 
pp. 54-62, 1930. 
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Consider again the set of pure states in the lower automaton given by equation 
(54). Let pj be the probability that the yth node is in the pure state \^ j {t^ . Then 
the probability density operator is defined by: 

5 /'(0=Zp>|^.(0){^.(0| 

J 

(56) 

and can be used to determine the probability of a particular quantum state. 
Differentiating equation (56) and substituting equation (40) gives the differential 
10 operator: 



ih—p = H'p-p'H 

(57) 



15 The operator on the right-hand side of equation (57) is termed the "commutator" and 
can be written as [i/, p]. Thus equation (57) becomes: 



20 



(58) 



Equation (58) characterizes the computation carried out by the upper automaton 1412 
in Figure 14. The results obtained by the upper automaton 1412 is the probability of 
the node being in the state | (r)^ . 

Given the density operator in equation (56) the ensemble average for the 
25 lattice can be determined by the trace of the matrix p C and is given by the 
expression: 
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{C) = trace{p'C) 
(59) 

where C can be any quantum mechanical operator. For example, the vibrational 
5 Hamiltonian characterizing the vibrational motion of the nodes is substituted for the 
operator C to give the ensemble average of the vibration states. The ensemble 
average calculated for an operator C is carried out by the transducer 305 in Figure 3. 

Notice that the state function 14^^. (/)^ is given as a chattering combination in 

equation (51) and can equivalently be written as a linear combination of the states in S 
10 as follows: 

l^.(0>=Z-.(0|^JO) 

(60) 

15 where the coefficients (/) are the chattering coefficients given by equation (54). 
The equivalence between the linear combination in equation (60) and the chattering 
combination in equation (51) is established in the following version of the Chattering 
Lemma: 

Let 5 = 11^^. (t))|^e//^,A^^ is finitej and let S be the set of chattering 

20 combinations of S. Let a real positive number s be given. There exists state 

fimctions |0^^ e 5 , defined for each tuple given by: 
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max 



f{|0>W>-|«>J%.(-))} 



<€ 



(62) 



for all . 

5 The Chattering Lemma says that every chattering combination of the form of 

equation (51) on a set S can be realized as a linear combination of elements from S 
with an arbitrarily small error. Under strong continuity assumptions of the state 
functions, the converse of the lemma is also true. By choosing the boundary 
conditions in the lattice appropriately, the continuity assumption is not limiting. 
10 Although the linear combination given in equation (60) can be used to 

represent the state function {t^ of the yth node, formulating the state function in 

the form of the chattering combination in equation (51) allows for formulating the 
sequence of excitation steps. 



15 3. The Input Selector Function 

The Input Selector Function 1418 in Figure 14 utilizes the information 
received from the current state 1410 of the node, the probability state 1414, models 
the interaction of a node in the lattice with the excitation field 1418 and with the 
electric fields generated by neighboring nodes 1420 to generate Hf. Figure 17 shows 

20 examples of the kinds of lattice node interactions. For the sake of simplicity, only 
nearest neighbor interactions are shown in Figure 1 7, however, the present invention 
can accommodate other node interactions. Oval 1702 identifies a dipole moment 
interaction resulting from a node 1704 interacting with a nearest neighboring node 
1706. Oval 1708 identifies a tripole moment interaction created by a node 1710 

25 interacting with two neighboring nodes 1712 and 1714. Oval 1716 identifies a 
quadrupole interaction between a node 1718 neighboring nodes 1720, 1722, 1724, and 
1726. 

In order to distinguish the excitation field from the electric fields generated by 
other nodes in the lattice, the excitation-field Hamilitonian Hf can be written as the 
30 sum of two terms as follows: 



I 
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where He is the interaction of the node with the excitation field, and Hi is the 
interaction of the node with other nodes in the lattice. The interaction Hamiltonian Hi 
for a given node can be written as: 



where for the yth node, H\^ is the interaction Hamiltonian of the node j with 
"neighboring" nodes i}, and is referred to as the "neighboring function" at the 
current state of the node | ^ . {t^ . Figure 18 illustrates the notation for a hypothetical 
interaction Hamiltonian Hi given by: 



Edges 1801-1803 identify the interaction Hamiltonians H\^ , H]' , and H]^ between 
the nodey and the neighboring nodes ly 1804, 2/ 1805, and 3y 1806, respectively. The 
neighborhood function represents the locus structure of interaction of the node j in 



The structure of interaction Hamiltonian Hi can be determined by constructing 
the Lagrangian associated with the lattice and applying the canonical quantization 
procedure discussed above in subsection I. D. The elements of carrying out this task 
are well known but the details depend on the physical characteristics of the lattice. 
The approach consists of defining potentials to characterize the lattice node 
interactions. For example, a dipole interaction as shown by oval 1702 in Figure 17 
may be characterized by a potential of the form: 



(64) 
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10 



(65) 



5 where d^_y (t) is a relaxation fiinction and nj is the distance between nodes / and J. 

For purposes of analysis and to determine a detailed formulation of the 
realization step described above in section I.F, assume that a finite set of primitive 
excitation Hamiltonians is given by: 



T = {H'^\i = l,...,N,} 



(66) 



where Ng is the number of excited states available to the node / The primitive 
excitation Hamiltonians given in equation (66) can be realized and implemented by 
15 chattering among the elements in the set rover the update time A. Chattering among 
the elements of Tis given by the chattering equation: 



20 



//>-|T(r)) 



(67) 



25 



where tiie sets are defined by equation (52) above and is a state 

function. The probabilistic resonance chattering equation (67) is central to the 
implementation of the quantum processor. The idea is to induce a probability 
distribution p on the states of the nodes in the lattice so that the realization criterion is 
satisfied. 
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Figures 1 9 and 20 show how a hypothetical quantum processor can be used to 
obtain an approximation to the function <P. Figures 19A-C illustrate a single 
hypothetical run of the excitation-field Hamiltonian Hf used to obtain a hypothetical 
spectrum from the quantum processor. In Figure 19 A, all of the nodes of the quantum 
processor 1902 begin in the ground state and can be characterized by the hardware 
Hamiltonian Hq. The detector 1904 coverts the incident radiation into an electrical 
current which is used to count the number of photons. In Figure 19B, the excitation- 
field /^jr impinged upon the quantum processor 1802 via the excitation generator 303 
in Figure. The excited nodes emit light radiation in the form of photons that strike the 
detector 1904, as indicated by radiating lines 1906, 1908, and 1910. Figure 19C 
shows a plot of the frequencies along the horizontal axis 1914 and the corresponding 
intensities along the vertical axis 1916 for the single application of the excitation-field 
Hamiltonian Hf, The radiating lines 1906, 1908, and 1910 are associated with nodes 
that emit radiation of the same frequency Vk and are included in determining the 
intensity value Xk, The resulting spectrum intensities {/tij/J,,...} shown in Figure 19C 
are stored in the coherent memory 306 in Figure 3. 

Figure 20 shows example hypothetical spectra and the average spectrum after 
n hypothetical runs of the excitation-field Hf. The n spectra indicated by spectrum 
2002, 2004, and 2006 are used to compute an average spectrum 2008. Approximately 
10^^ runs are used to determine the average spectrum 2008 where each run takes about 
.2 nanoseconds to complete. Rather than store the spectrum for each run, a rurming 
average of the intensities, referred to as [aJ,^,...}, is stored in the coherent memory 

306 in Figure 3. The average spectrum can be used to approximate the function 
0(jc) in equation (42) by the following: 

(68) 



D. The Read-Out Period 



r 
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The read-out period, referred to as is the amount of time the transducer 
devotes to reading the light radiation emitted from the quantum processor in step 410 
in Figure 4. If the transducer collects singles from the quantum processor for too 
long, there is a danger of de coherence due to emissions caused by nodal interaction 
5 characterized by the interaction Hamiltonian ///. For example, a computation carried 
out in the lattice of the quantum processor begins with a propagating probabilistic 
wave train as described above in relation to Figure 13. Figure 21 is an illustration of a 
forward propagating wave train 2102 and a backward propagating wave train 2104, 
After the initial excitation-field i:^ impinges upon the lattice, the forward propagating 

10 wave train 2102 sweeps across the lattice of nodes, where some of the energy stored 
in the excitation field is converted into nodal excitation energy. The nodes respond 
after excitation by relaxing to lower energy states by releasing light radiation hv2\0%. 
The excitation field sloshes back and forth in the form of forward and backward 
moving wave trains 2102 and 2104, respectively, until the entire excitation field Hf 

15 has been consumed and converted to light energy via nodal excitation. During this 
time period radiation continues to be emitted as a result of nodal interaction 
characterized by the interaction Hamiltonian Hi of different energies A v'21 10, which 
can distort the intensity-versus-frequency spectrum. Therefore, in order to limit the 
amount of light radiation emitted as a result of nodal interaction from contributing 

20 significantly to the intensity-versus-frequency spectrum, the readout period ^ for the 
transducer is considerably shorter than the amount time it takes for the excitation field 
to dissipate. 

Figure 22 is a control-flow diagram depicting determination of the transducer 
read-out period ^ representing one of many possible embodiments of the present 

25 invention. In step 2202, the quantum processor lattice is at an initial probabilistic 
state characterized by the hardware Hamiltonian Hq. In step 2204, the programmed 
excitation-field //^determined by the ftinction input 302 and realized by the excitation 
generator 303 in Figure 3 impinges upon the quantum processor lattice of nodes. In 
step 2206, after a read-out period^' has elapsed, the transducer 305 in Figure 3 

30 transmits the spectrum intensities and corresponding frequencies to the coherent 
memory 306 in Figure 3. In step 2208, if the resulting equation (68) does not 
approximate the encoding ftinction within parameters set by the operator, then in step 
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2210, the read-out period ^ is extended to a longer read-out period ^/ and the 
computation is started again from the initial probabilistic state. In step 2208, if the 
resulting average spectrum intensities are satisfactory, then the 

computation is complete and an approximation to the encoding function has been 
5 determined. 

The read-out period ^ cannot be extended arbitrarily because thermal 
relaxation mechanisms inherent in a quantum mechanical system will eventually 
induce de coherence. In order to demonstrate this, consider the discrete spectrum 
[E.] of the hardware Hamiltonian Ho given by the eigenvalue equation: 

10 

//o|4>,(0> = £,|%(0> 

(69) 

where |^/(^)) are the corresponding eigenfunctions. Using equation (58), the 
15 transition probability matrix elements p.j satisfies the expression: 

(70) 

20 where the matrix elements are given by: 

0 

(71) 

25 where *(^,0 complex conjugate of the function (;c,/) . If the system is 

in the state at time / equal to zero, then the presence of the relaxation 

Hamiltonian Hi given by equation (54) causes the corresponding probability density 
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matrix elements to decay exponentially with time. See Fields and Nodes: Field 
Theory and Dispersion Relations, K. Nishijima, Benjamin, Inc., N.Y., 1969. For 
small values of t, the diagonal matrix elements Pjj are the largest terms in the matrix 
representation of the probability density operator. Assuming no excitation and using 
5 equation (57), the term satisfies the following equation: 



15 



(72) 



10 where 



dt 
(73) 

for / not equal to j. The solution of equations (72) and (73) is given by: 



/4 

(74) 



where ^ is the relaxation time that depends on the physical characteristics of the 
20 quantum processor lattice. The relaxation time is the time required for the wave train 
to achieve de coherence. Therefore, the effect of the inter-node excitation 
Hamiltonian Hi is to randomize the state transitions in each node corrupting the 
computation described above . Therefore, the read-out period ^ must be chosen so 
that this effect is acceptable. In other words, the read-out period is selected so that 



(75) 
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The upper bound ^ is used in the design specification of the quantum processor and in 
the computability analysis. 

Although the present invention has been presented in terms of a particular 
embodiment, it is not intended that the invention be limited to this embodiment. 
Modifications within the spirit of the invention will be apparent to those skilled in the 
art. 

The foregoing description, for purposes of explanation, used specific 
nomenclature to provide a thorough understanding of the invention. However, it will 
be apparent to one skilled in the art that the specific details are not required in order to 
practice the invention. The foregoing descriptions of specific embodiments of the 
present invention are presented for purposes of illustration and description; they are 
not intended to be exhaustive or to limit the invention to the precise forms disclosed, 
obviously many modifications and variations are possible in view of the above 
techniques. The embodiments were chosen and described in order to best explain the 
principles of the invention and its practical applications and to thereby enable others 
skilled in the art to best utilize the invention and various embodiments with various 
modifications as are suited to the particular use contemplated. It is intended that the 
scope of the invention be defined by the following claims and their equivalents: 



